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Abstract 

We propose a simple pedagogical way of introducing the Euler-MacLaurin summation formula in 
an undergraduate course on statistical mechanics. We put forward two alternative routes: the first 
one is the simplest and yields the first two terms of the expansion. The second one is somewhat 
more elaborate and takes into account all the correction terms. We apply both to the calculation 
of the simplest one-particle canonical partition functions for the translational, vibrational and 
rotational degrees of freedom. 
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I. INTRODUCTION 


In statistical mechanics the thermodynamic fnnctions are given in terms of the logarithm 
of the partition function and its derivatives with respect to volume, temperature, etc. The 
mathematical expression of the canonical partition function is an inhnite sum over all the 
states of the system or, under some simplifying assumptions, the states of the individual 
particles.- Some of these inhnite sums are commonly calculated approximately by means of 
the Euler-MacLaurin summation formulai^ and most textbooks simply show how to apply it 
to the cases of interest.- A rigorous derivation of the summation formula may be rather too 
demanding for inexperienced students and one avoids it in introductory courses on statistical 
thermodynamics or statistical mechanics. 

There is a relatively simple way of deriving the summation formula by means of oper¬ 
ator methods^ but it also requires some kind of mathematical expertise that students of 
introductory courses may not posses. However, the students may feel more conhdent and 
confortable if they are shown how to carry out such calculations by means of mathematical 
methods that they have already learned in a hrst course on mathematical analysis. 

In an introductory undergraduate course we show how to obtain the hrst two terms of 
the EuIer-MacLaurin summation formula in an extremely simple way that only requires the 
students to be familiar with the Taylor expansion. The aim of this paper is to put forward 
this approach that we deem suitable for pedagogical purposes. In section [TTl we outline the 
method and apply it to some simple examples in section HTTl If we decide to show the students 
how to obtain terms of higher order this simple approach results to be rather cumbersome. 
For that reason, in section nYi we show how to derive the full summation formula by a 
relatively minor change of strategy that requires some additional mathematical skills. The 
reader may choose one or another approach depending on the level of the course. Finally, 
in section |V] we summarize the results and draw conclusions. 


II. THE SIMPLEST APPROXIMATION 

The simplest canonical partition functions in statistical mechanics can be expressed as 
sums of the form 

CO 

S = J2f(n), (1) 

n=Q 
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where /(n) should tend to zero sufficiently fast when n —?• cx) because, otherwise, the sum 
does not converge. In most cases one cannot obtain this sum in closed form and therefore 
resorts to some kind of approximation valid under certain conditions; for example, sufficiently 
high temperature. In order to derive such an approach we dehne 

OO 

S{x) = '^f{n + x), (2) 

that satishes 

^(0) = 5 

lim S{x) = 0, (3) 

X—)-oo 

and 

S{x)-Six + l) = fix). (4) 

Substituting the Taylor expansion S{x + h) = S{x) + S\x)h + ^S"{x)h‘^ + ... for h = 1 
into (jl]) 

we obtain, after some rearrangement, 

- S'(x) = /(x) + is"(x) + .... (5) 

If we integrate this expression between x and oo and take into account that S{x) and all its 
derivatives vanish at the upper limit we have 

1 

S{x) = j -S\x) + (6) 

that can be solved iteratively. In the hrst step we omit the derivatives of S(x) so that this 
function is approximately given by S(x) ~ If we substitute this approximate 

result into the right-hand side of (E]) we have 

f{t)dt+^f{x). (7) 

Thus, for X = 0 we obtain a simple approximation to the sum ([1]): 

S(0) = S« r/(i)o!i+l/(0). (8) 

The simple method just outlined is suitable for introductory courses because it only 
requires basic knowlegde in mathematical analysis. This procedure is not suitable for the 
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systematic calculation of the corrections of higher order to the basic formula (j8]) because it 
soon becomes rather cumbersome. However, this result is suitable for most porposes in an 
introductory course on statistical mechanics because several partitions functions for atomic 
and simple molecular systems can be easily derived from it.^ 


III. SOME SIMPLE EXAMPLES 


Under some simplifying assumptions that we will not discuss in this paper statistical 
mechanics tell us how to express the thermodynamic functions of atomic and molecular 
systems in terms of one-particle partition functions of the form^ 






( 9 ) 


where V is the volume of the container, T the absolute temperature, the n-th energy level 
(assumed to be ^fri-fold degenerate) of the particle and ks the Boltzman constant. 

Let us hrst consider a particle of mass m in a one-dimensional box of length L with 
impenetrable walls which is the starting point for the calculation of the thermodynamic 
properties of an ideal monoatomic gas.- If the particle has no internal structure its spectrum 
is only given by the translational degree of freedom: 

( 10 ) 




where h is the Planck constant h divided 27r and Qn = 1. Upon dehning a = h^/( SmkBTL'^) 
the partition function qt reads 

qt = St-I 


and the approximation ([8]) yields 


St = 

(11) 

n=Q 


St - ^ + ^ 

* 2^^ 2' 

(12) 


This approximate expression is accurate for sufficiently small values of a or, equivalently, 
when the thermal de Broglie wavelength A = h/y/2'KmkBT is much smaller than the box 
length L.l The addition of corrections of higher order commonly improves the result of the 
approximate partition function, but in this case all of them vanish.- Figure [1] shows that 
this approximate expression becomes increasingly more accurate as a decreases. 
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The spectrum of a one-dimensional harmonic oscillator is given by 

el= + n = 0,1,..., (13) 

where v is the frequency of the oscillation and gn = The partition function for this degree 
of freedom is given by 


n — 

oo 

(14) 

n=0 


where a = hu/{kBT). In this case the sum is a geometric series that can be calculated 
exactly and therefore we can derive its small-a expansion in a straightforward way: 

C ^ ^ I ^ I “ “^1 

“ 1 - e-" ~ a~^ 2~^ 12~ "" 



If we define the vibrational temperature 6y = hv/hs then a = 9^/T and the Euler- 
MacLaurin formula is a good approximation when 9y <C T. The approximate expression ([8]) 
yields the hrst two terms of this series exactly. 

The last example is the partition function for a rigid rotor with moment of inertia I. In 
this case, every energy level 

= = ( 16 ) 


is (2J -|- l)-fold degenerate and the partition function for the rotational degree of freedom 
is given by 

OO 

qr = Sr = (17) 

J=0 

where a = / {2IkBT). The approximate expression ([8]) yields 


Sr 



(18) 


that is accurate enough if a is sufficently small as shown in hgure [2l It is custommary to 
dehne the rotational temperature 9r = / {2IkB)r so that a = 9r/T 1 when 9r -C T. It 

is worth mentioning that the a-independent term in equation (lT8ll is not exact as shown in 
the following section. 
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IV. SYSTEMATIC APPROACH 


If we decide to show the students how to obtain expression of higher order it is convenient 
to proceed in a different way. If we write the Taylor expansion discussed in section [ITl as 

OO 

sn+i) = ^-iS«(n, (19) 

J-o 

where S^^\x) is the j-th derivative of S{x) with respect to x, then equation ([6]) reads 

OO ^ 

S(n = f (i) - XJ 

where F{x) = f{t) dt. 

Instead of trying to solve equation fl20l) iteratively we propose a solution of the form 

OO 

= ( 21 ) 

k=0 

where F^^^x) = x). In order to obtain the coefficients we substitute fl2l|) into 

(|20|) and compare the coefficients of F^^\x) in the left- and right-hand sides; the result is 

n 

E ^n—j „ 

(j + 1)1’ ^ - 1)2, , 

ao = 1. (22) 


We thus have 

OO 

^(x) = F{x) - ak+if^’'\x), 

k=0 

and 

POO ^ 

S= / /(t)dt-5^afc+i/(^)(0). 

k=0 

Straightforward inspection of the hrst coefficients a„ 


Oi 


as 


~2’ ~ 12’ — 0) — 0, 06 

1 _n _ ^ 

~ 1209600’ 1 «io - 47900160’ 


1 

30^’ ^ 


(23) 


(24) 


(25) 


suggests that 02^+1 vanish for all j > 0. In order to prove this conjecture we consider the 
function 


u{x) = 


X 


- 1 ' 


(26) 
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that satisfies 

u{x) — u{—x) = —X. (27) 

If we substitute the Taylor expansion 

OO 

u(x) = ^^UjX^, (28) 

j=0 

into equation fl^ we conclude that U2j+i = 0 for all j > 0. If we now expand (e®—l)n(a;) = x 
in a Taylor series about x = 0 and compare the coefficients of x” in the left- and right-hand 
sides of the resulting equation we obtain a recurrence relation for the coefficients Un that is 
identical with equation for the coefficients a„. We thus conclude that for all 

n. Therefore, equation fl2T|) reduces to 

^OO -j 

S= / /(i)*+ 2/(0)-5^a2t/'“-‘>(0). (29) 

k=l 

It is worth noting that the coefficients aj are related to the Brillouin numbers Bj^ in the 
following way: aj = Bj/jl. 

The operator method leads to this result in a more straightforward way and is also more 
convenient for the discussion of the radius of convergence of the series.- However, we do not 
discuss it here because it requires the introduction of functions of operators that may not 
be suitable for an undergraduate course. 

Let us apply this more accurate approach to the examples discussed in section IIIII As 
pointed out in that section nothing can be done with the particle in a box because all the 
corrections of higher order vanish. This surprising result can be explained very easily and 
applies to all the sums in which f(—x) = f(x) because x) = —and 

y(2fc-i)('g^ = 0. A more rigorous analysis of such problems, based on the Poisson summation 
formula, can be found eslewhere.- 

The harmonic oscillator is a suitable simple example for testing the summation formula 
fl29|) because we can calculate the exact expansion as shown in equation flT^ . In this case 
every new term added to the summation formula 0291) yields one more term of the small-a 
series (lT5l) as one may easily verify. 

The application of the summation formula fl2^ to the partition function for the rigid 
rotor should be carried out with care. The reason is that = a^Pk+i(a), where 

Pk+i(a) is a polynomial function of a of degree k + 1 and Pk+i(0) ^ 0 . Therefore, the 


7 


summation formula that includes all the terms through only yields the small-a 

series correctly through degree k. For this reason the simple expression (|H]) does not yield the 
correct a-independent term that receives contributions from /(O) and f'{0)- The expansion 
accurate to is easily shown to be 


Sr 


1 ^ 1 ^ a ^ ^ ^ ^ 

o^3^15’''315^315^ 3465 


(30) 


if we add all the terms trough /^®^(0) in the summation formula fl2^ as argued above. Since 
this series does not converge^ one should truncate it before the terms start to increase. 


V. CONCLUSIONS 

As stated in the introduction we show the students the simple method developed in sec¬ 
tion [ITl because it is sufficient for the purposes of our course. However, if a motivated student 
wants to learn how to obtain the corrections of greater order that appear in some of the 
available texbooks on statistical mechanics then one can suggest him or her to try the sys¬ 
tematic approach of section [I^ or even the operator method.- In our opinion the present way 
of deriving the Euler-MacLaurin summation formula for statistical mechanics applications 
is more convenient for pedagogical purposes than the traditional one that appears in most 
textbooks on mathematics or numerical analysis.- Such traditional approaches are certainly 
more rigorous and general but require the students to be more experienced in mathematics. 

In passing, it is worth mentioning that the summation formula given by equations (122|1 
and fl29l) is suitable for motivating the students to resort to a computer algebra system. In 
this way they bypass the tedious algebraic manipulation of the equations that is necessary 
for the calculation of contributions of large order like those in equation fl30|) and practise 
programming in any of such useful languages. 
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FIG. 1: St calculated by means of the sum (11111 (points) and the simple Euler-MacLaurin expression 
(fT^ (solid line) 
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FIG. 2: Sr calculated by means of the sum (fTTll (points) and the simple Euler-MacLaurin expression 
(fT8|) (solid line) 
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